WELL-POSEDNESS THEORY 
FOR GEOMETRY COMPATIBLE HYPERBOLIC 
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Abstract. Motivated by many applications (geophysical flows, 
general relativity) , we attempt to set the foundations for a study of 
entropy solutions to nonlinear hyperbolic conservation laws posed 
on a (Riemannian or Lorentzian) manifold. The flux of the conser- 
vation laws is viewed as a vector-field on the manifold and depends 
on the unknown function as a parameter. We introduce notions 
of entropy solutions in the class of bounded measurable functions 
and in the class of measure- valued mappings. We establish the 
well-posedness theory for conservation laws on a manifold, by gen- 
eralizing both Kruzkov's and DiPerna's theories originally devel- 
oped in the Euclidian setting. The class of geometry-compatible 
(as we call it) conservation laws is singled out as an important case 
of interest, which leads to robust L p estimates independent of the 
geometry of the manifold. On the other hand, general conservation 
laws solely enjoy the L 1 contraction property and leads to a unique 
contractive semi-group of entropy solutions. Our framework allows 
us to construct entropy solutions on a manifold via the vanishing 
diffusion method or the finite volume method. 



1. Introduction 

The theoretical work on discontinuous solutions to nonlinear hy- 
perbolic conservation laws has been restricted so far to problems set 
in the Euclidian space MJ 1 . Motivated by numerous applications, to 
geophysical fluid flows (shallow-water equations on the surface on the 
Earth) and general relativity (Euler-Einstein equations describing neu- 
tron stars and black holes) in particular, we attempt in the present 
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paper to set the foundations for a study of weak solutions (including 
shock waves) to hyperbolic conservation laws posed on a Riemannian or 
Lorentzian manifold. Recall that Kruzkov's theory [6] deals with equa- 
tions posed in the Euclidian space M n and provides existence, unique- 
ness, and stability of L 1 D L°° entropy solutions to nonlinear hyperbolic 
balance laws in several space dimensions. In addition, in the case of 
conservation laws with "constant flux" , depending solely on the conser- 
vative variable but not on the time and space variables, solutions are 
known to satisfy the maximum principle as well as the L 1 contraction 
and total variation diminishing (TVD) properties. 

The above stability features play a central role in the theory and 
numerical analysis of hyperbolic conservation laws. The purpose of 
this paper is to introduce a suitable generalization of Kruzkov's theory 
when solutions are defined on a manifold rather than on R™. An investi- 
gation of the interplay between the geometry of manifolds with limited 
regularity and the analysis of hyperbolic partial differential equations 
(for non-smooth functions) appears to be particularly challenging. We 
hope that the framework initiated here will lead to further works on 
this important subject. 

Let (WL n ,g) be an oriented, compact, n-dimensional Riemannian 
manifold. One of our tasks will be to define a suitable class of con- 
servation laws posed on M n to which the stability properties of the 
Euclidian space could possibly extend. It should be observed that, on 
a manifold, there is no concept of "constant flux" -in the sense that 
the flux would be "independent" of the space (and time) variable, since 
as we will see below the flux is a section of the tangent bundle TM n . 
Indeed, the flux at a point x G M n belongs to the tangent space T^M™ 
to the manifold at that point. The Euclidian space is special in that it 
is possible to choose the flux f(u) at an arbitrary point x G W 1 and to 
parallel-transport it to the whole space M n ; the vector field generated 
in this fashion remains constant in x and, therefore, smooth. By con- 
trast, if a tangent vector is selected at one point of the manifold one can 
not, in general, associate with it by parallel-transport a unique smooth 
vector field on the manifold; this is a consequence of the curvature of 
M n . As an example, we recall that any vector field on S 2 (the unit 
sphere in R 3 , an important example arising in geophysical flows and 
motivating our study) must have two critical points, at least, or else 
must be discontinuous. 

An outline of this paper is as follows. In Section [3j we introduce 
the class of geometry-compatible conservation laws on M" based on 
divergence-free flux, which we single out as a class of particular inter- 
est. We then define the notions of entropy pair and entropy solution 
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in the function space L°°(IR + x M"). In Section HJ we show the exis- 
tence of entropy solutions via the so-called vanishing diffusion method, 
and establish that solutions to geometry-compatible conservation laws 
enjoy LP stability properties (for all p) that are independent of the 
geometry of the manifold. Next, Section [5] is devoted to the general 
well-posedness theory on a Riemannian manifold, and covers both L°° 
solutions (following Kruzkov |6j) and measure- valued solutions (follow- 
ing DiPerna [3]). 

In turn, we obtain a versatile framework allowing us to establish 
the convergence toward the unique entropy solution of any sequence of 
approximate solutions satisfying all of the entropy inequalities (possibly 
up to some small error). Our geometry-independent bounds should be 
useful for the numerical analysis of stable and robust, shock-capturing 
schemes adapted to hyperbolic equations posed on manifolds. In this 
direction, the follow-up paper [1] will cover the convergence of finite 
volume methods based on non-Cartesian meshes, and the derivation of 
total variation diminishing (TVD) estimates. 

The remaining of the paper contains several generalizations of in- 
terest. In Section [HI we show that the well-posedness theory can be 
developed in the function space L l for general conservation laws that 
need not be geometry-compatible; we prove that the L 1 contraction 
property hold at this level of generality, but that for p > 1 the LP 
stability properties may be violated. Finally, in Section [7] we discuss 
the well-posedness theory for conservation laws posed on a Lorentzian 
manifold, including for instance the Schwarzschild spacetime. 



2. Preliminaries 

In this preliminary section, we briefly review some key concepts and 
results from Kruzkov's theory in the Euclidian space. Consider the 
hyperbolic conservation law 

n 

d t u + ^2dj(f j (u,-)) = 0, u = u{t,x) eR,t>0, x eR n , (2.1) 

3=1 

where dj := d/dx^ and / : R x MP — > M. n is a given smooth mapping, 
refered to as the flux of (12.11) . Observe that in (12.11) (as well as in (12.31) 
below) the divergence operator acts on the x-coordinate and takes into 
account the dependence in x = (x J ) of both the solution u = u{t,x) 
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and the flux / = (f J (u,x)), that is 

n 

3=1 

n 

= U d uf){u{t,x),x)dju{t,x) + (djf J )(u(t,x),x)y 

3=1 

As is well-known [21 EJ El CO El EG] , solutions of (12.11) are generally dis- 
continuous and must be understood in the weak sense of distributions. 
In addition, for the sake of uniqueness one should restrict attention 
to entropy solutions characterized by an infinite family of inequalities. 
For simplicity, assume first that the flux is divergence- free, that is 

n 

i d jf j ) fa x) = 0, ueR,xe R n . (2.2) 

3=1 

(This assumption is motivated by Lemma [2.11 below.) When the con- 
dition (12.21) is imposed, the entropy inequalities associated with the 
equation (12. ip take the form 

n 

d t U(u) + Y,d 3 {F^u,-))<0, 

3=1 

U" > 0, 

/u 
d u U(u) d u f J (u, x) du', x G M n , u e R, j = 1, . . . , n. 

(2.3) 

Relying on Kruzkov's theory [6] and the arguments given later in this 
paper, one can check that, given any divergence-free flux and initial 
data m e L^M 71 ) n L°°(R n ), the initial value problem (ETC]) . 

u{0,x) =u {x), xeR n , (2.4) 

admits a unique entropy solution u £ L°° (R + , L 1 (M ri )nL 00 (R' 1 )) which, 
moreover, satisfies the LP stability properties 

IK*)||lp(r«) < ||«o||lp(r»), t > 0, p e [l,oo], (2.5) 

together with the L l contraction property: for any two entropy solu- 
tions u, v 

\\v(t) - u(t) || L i (R n) < \\v(t') - u(t') \\ L x m , 0<t'<t. (2.6) 

Furthermore, one has the following regularity property: If the initial 
data has bounded total variation, that is Uq G BV(R n ), then the solu- 
tion has bounded variation for all times, that is u G L°°(lR + , BV(R n )). 
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In the Euclidian setting, the condition ( 12.2ft arises in the following 
way. 

Lemma 2.1. A necessary and sufficient condition for every smooth 
solution of (12.11) to satisfy the infinite list of additional conservation 
laws (12. 3p is that the divergence-free condition (12. 2p holds. 

Proof. This follows readily by multiplying the conservation law (12.11) 
by the derivative d u ll(u) of an arbitrary function U(u), using the chain 
rule, and taking advantage of the fact that the function U is arbitrary. 
Namely, we have 

n n 

dj(F j (u, ■)) = J>„^')(«, ■) d 3 u + (d^)(u, ■) 

3=1 3=1 
n 

= {d u Ud u f)(u, ■) d jU + (djF^u, ■), 

3=1 

so that 

n 

d t U(u) + dj(Fi(u, ■)) = ■) - d u U{u) (d^u, ■). 

3=1 

Now, by imposing that the right-hand side of the above identity van- 
ishes for all entropy U and by differentiating in u the corresponding 
relation while using the definition of in ( 12. 31) . we obtain 

d uu U{u) (d j f j )(u,-) = 0. 
Since U is arbitrary, the desired claim follows. □ 

Observe that Kruzkov's theory in M. n also covers flux that need not 
satisfy our condition (I2.2p . and actually applies to balance laws of the 
general form 

n 

d t u + J2 d j{f j ( u r))=h{u,-), u = u{t,x)eR, (2.7) 

3=1 

where / = (ft(u,t,x)) and h = h{u,t,x) are given smooth mappings. 
Provided h is Lipschitz continuous in the u- variable (uniformly in t, x), 
the initial value problem associated with ( 12.71) admits a globally de- 
fined, unique, entropy solution u G L^ C (IR + x W 1 ). Note, however, that 
the general equation (12.71) does not enjoy the LP stability, L l contrac- 
tion and TVD properties. Most of the literature on scalar conservation 
laws is restricted to the case of a constant flux / = f(u) which, clearly, 
arises as a special case of (12.21) . As pointed out earlier, there is no such 
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concept as a "constant flux" in the context of manifolds, and a suitable 
generalization of f)2.2p will arise in our analysis. 

3. Geometry-compatible conservation laws 

Background and notation. Throughout this paper, (NL n ,g) is a 
compact, oriented, n-dimensional Riemannian manifold. As usual, the 
tangent space at a point x G M n is denoted by T x .M n and the tan- 
gent bundle by TM™ := [j xeM n T x ~M n , while the cotangent bundle is 
denoted by T*M n = T*M n . The metric structure is determined by a 
positive-definite, 2-covariant tensor field g, that is, at each x G M n , g x 
is a inner product on T X M™. 

In local coordinates x = (x^), the derivations dj := ■A- form a basis 
of the tangent space T x Nl n , while the differential forms dx^ determine a 
basis of the cotangent space T*M n . Here and below, we use Einstein's 
summation convention on repeated indices so, for instance, in local 
coordinates 

g = g,^ dx l dx\ g^ = g(d h dj). 

The notation < X x , Y x > gx '-= g x (X x , Y x ) may also be used for the inner 
product of X x , Y x G T X M™. 

We denote by (g^) the inverse of the positive definite matrix (gij). 
The metric tensor will be used to raise and lower indices, so that to 
each vector X = one associates the covector (Xj) via 

Xj := g^ X\ 

Recall that the differential df of a function / : M n — » R is the section 
of the cotangent bundle T*M n , that is df x G T*M n , defined by 

df x (X x )=X x (f), X x eT x M n . 

We denote by dV g the volume form on M ra , which in local coordinates 
reads 

dVg — \f\g~\dx 1 ■ ■ ■ dx n , 

with \g\ := det(^y) > 0. The gradient of a function h is obtained 
from dh by using the isomorphism between TM" and T*M n via the 
Riemannian scalar product. In local coordinates, the gradient of a 
function h is the vector field 

f ^h={Vh) = ± f :=g«8 i h±. 

To the metric tensor g one associates a covariant derivative operator 
V, characterized by the condition Vg = 0. In particular, the covariant 
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derivative of a vector field X = X- 7 is the (1, l)-tensor field whose 
coefficients in local coordinates are 

V,,V : djX* • l^.V'. i,j = l,...,n, 

where the Christoffel's symbols are determined from the metric tensor 
by 

The divergence of the vector field X is the function 

<Hv,.V := VjXi = djX* + 1% X k . (3.1) 
Recall the duality relation 

/ g(grad g h,X)dV g = - I hdiv g XdV g , (3.2) 

which is valid for any smooth function h and vector field X, at least. 
Interestingly enough, it follows from (13.21) that, in local coordinates, 

div 9 X = (yfigi^djiVWlX*), (3-3) 
which shows that the divergence operator on vector Gelds depends upon 

\g\, onl y- 

A class of conservation laws. The observation made earlier in Lem- 
ma l2.1l in the Euclidian setting provides us the motivation for the defini- 
tion that we now introduce. As was already pointed out in the introduc- 
tion, although Kruzkov's theory applies to general balance laws (12. 7p . 
the subclass defined by (12. ip and the condition (12.21) leads to many 
properties (maximum principle, L 1 contraction, L p stability) which are 
very desirable features that we attempt to guarantee on the manifold 
M n . At this stage of the discussion, we restrict attention to smooth 
functions satisfying in a classical sense the partial differential equations 
under consideration. 

Definition 3.1. 1. A flux on the manifold M n is a vector field f = 
f x (u) depending upon the parameter u (the dependence in both variables 
being smooth). 

2. The conservation law associated with the flux / on M" is 

d t u + div g (f(u)) = 0, (3.4) 

where the unknown is the scalar-valued function u = u(t, x) defined for 
t > and x G M n and the divergence operator is applied, for each fixed 
time t, to the vector field x \— > f x (u(t,x)) G T x JsA n . 

3. A flux is called geometry-compatible if it satisfies the condition 

div g f x (u) =0, u G R, x G M n . (3.5) 
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We will also refer to ( I3.4p -( l3~5j) as a geometry-compatible conserva- 
tion law. Let us emphasize that (13.4)) is a geometric partial differential 
equation which depends on the geometry of the manifold, only, and 
is independent of any particular system of local coordinates on M. 
In particular, all estimates derived for solutions of (13.41) must take a 
coordinate-independent form, while, in the proofs, it will often be con- 
venient to introduce a particular chart. 

Observe that, by the divergence theorem on Riemannian manifolds 
([llj. for instance), all sufficiently smooth solutions of (13.41) satisfy the 
balance law 



for every smooth, n-dimensional sub-manifold S C M n with boundary 
dS. Here, v denotes the outward unit normal to S, and Vgs is the 
volume form induced on the boundary by restricting the metric g to 
the tangent space T(dS). 

In local coordinates, the flux components f x (u) depend upon u, x 
and we will denote by d u f ] x {u) and dkfi(u) its u- and ^-derivatives, 
respectively. Before we proceed, observe that: 

Lemma 3.2. Let f = f x {u) be a geometry- compatible flux on M n . In 
local coordinates the conservation law (13.41) takes the (nonconservative) 
form 



Proof. It follows from the definition of the divergence operator (13.11) 
that the conservation law (13.41) can be written as 



spect to the local coordinate x of the composite map x h- > f x (u(t,x)), 
that is 




d t u(t, x) + (d u f 3 x ) (u(t, x)) dju(t, x) = 0. 



(3.6) 




'Wi ("(/..<'))) = (d u fi)(u(t,x))d j u(t,x) + (d j fi)(u(t,x)). 



On the other hand, the condition in Definition 13.11 yields 



{Ojfi) {». x) + l%(x) f k x {u) = 0, u e R, x e M n . 



Writing this identity with u = u(t, x) and combining it with the above 
observations lead us to (13.6)) . □ 
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Weak solutions. We will be interested in measurable and bounded 
functions u G L°°(R + x M") satisfying (13.41) in the sense of distribu- 
tions and assuming a prescribed initial condition w -L°°(M n ) at the 
time t = 0: 

u(0,x)=u (x), xeM n . (3.7) 

For the sake of uniqueness, we need an analogue of the entropy in- 
equalities (12.31) proposed by Kruzkov in the Euclidian space R n . These 
inequalities are derived as follows. 

First of all, consider smooth solutions u = u(t,x) of (13.41) and mul- 
tiply the conservation law by the derivative d u ll(u) of an arbitrary 
function U. We obtain 

d t U (u(t, x)) + d u U (u(t, x)) dw g (f x (u(t, x))) = 0, 

which suggests to search for a vector field F = F x {u) such that 

d u U(u(t,x))div g (f x (u(t,x))) = divg (F x (u(t,x))). 

This relation should hold for all functions u = u(t,x), and is equiva- 
lent to impose that the flux components f ] x {u) satisfy the two partial 
differential equations 

d u U{u) d u fl(u) = d u F*(u), l<j<n, 

8 u U{u) (djfM + Y{ 3 {x) />;:(„)) = djFi(u) + I%{x) F';:(a). 

for all u G R and x G M™. The first equation implies that F x is given 
by 

F x (u):=J d w U(w)d w f x (w)dw eT x M". (3.9) 

Differentiating the second equation with respect to u and using (13.91) 
we find 

d uu U(u) (<) h fi(a) ■ Vip-)/^,,)) 0. 

Since U is arbitrary, a necessary and sufficient condition for (13.81) to 
hold is that f x {u) satisfies (13. 5p . Note in passing that necessarily F x (u) 
satisfies also this condition. 

This derivation provides us with an analogue of the definition (12.31) 
of the Euclidian case. In turn, smooth solutions of any geometry- 
compatible conservation law (13.41) satisfy the additional conservation 
laws 

d t U(u(t, x)) + div g {F x (u(t, x))) = 0, (3.10) 

where U is arbitrary and F is given by (13.91) . 

We are now in position to introduce the notion of entropy solution 
which generalizes Kruzkov's notion to the case of manifolds. We require 
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the equalities (13.101) to hold as inequalities only, as this arises naturally 
by the vanishing diffusion method. (See Section S]). 

Definition 3.3. Let f = f x (u) be a geometry- compatible flux on the 
Riemannian manifold (M,g). 

1. A convex entropy /entropy-flux pair is a pair (U,F) where U : 
R — > R is a convex function, and F = F x (u) is the vector field depend- 
ing on the parameter u defined by 



F x {u) := / d u ,U{u') d u ,f x {v!) du', u G R, x G M n . (3.11) 



2. Given a function u Q G L°°(M n ), afunctionu in L°° (R + , L°°(M n )) 
is called an entropy solution to the initial value problem (I3.4l) - (l3.7p if 
the following entropy inequalities hold 



for every convex entropy/ entropy flux pair (U, F) and all smooth func- 
tions 9 = 9(t,x) > compactly supported in [0, oo) x M n . 

4. Vanishing diffusion method on manifolds 

In the present section, we construct solutions of the conservation law 
(13.41) via the vanishing diffusion method. We establish the existence of 
smooth solutions to a regularized version of (13. Ah taking into account a 
small diffusion term, and then establish that these solutions converge to 
an entropy solution in the sense of Definition 13.31 as the diffusion tends 
to zero. For convenience, we provide a proof that is based on a uniform 
total variation bound, and refer the reader to the forthcoming section 
for an alternative approach based on the concept of measure-valued 
solution. 

Recall that M n endowed with the metric tensor g is a compact, ori- 
ented, n-dimensional Riemannian manifold. We denote by i/(M n ; dV g ) 
the usual Lebesgue spaces on the manifold (M™, g). The Sobolev space 
H 1 (M n ;dV g ) is the space of all functions h G L 2 (M n ; dV g ) such that 
grad fl h G L 2 (M"; dV g ) — where the gradient is defined in the distribu- 
tional sense, via the formula (I3.2p . The spaces if fc (M"; dV g ) for k > 2 
are defined similarly. 





U(u{t, x)) d t 9(t, x) + g x (F x (u(t, x)), grad g 9{t, x)) ) dV g (x)dt 




(3.12) 
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The total variation of a function h G L 1 (M n ; dV g ) is denned by 
TV(h) := sup / hdivgtp dVg, 

IMIi°°(M»)<l J m™ 

where ^ describes all C 1 vector fields on the manifold. We denote 
by BV(M n ; dV g ) C L l {M n ] dV g ) the space of all functions h with finite 
total variation, endowed with the norm ||/i||Li(M™;(iv )+TV(h). For ma- 
terial on BV functions we refer to [HH2]- m particular, it is well-known 
that the inclusion map of BV(Wl n ; dV g ) in L 1 (M n ; dV g ) is compact. 
It is clear that, for all smooth functions h : M" — > M, 

TV(h)= [ \&ad g h\ g dV g , 

where | ■ \ g denotes the Riemannian norm associated with g. Observe 
also that, by taking a partition of unity determining a finite covering 
of M™ by coordinate patches we can easily "localize" the concept of 
total variation on the manifold, as follows. 

Lemma 4.1. There exist finitely many (smooth) vector fields JW,. . . , 
on M n ; a constant Co > 1 (depending only on the manifold), and 
a chart covering the manifold such that: 

(ii) each vector field X® is supported in a coordinate patch and 

Span{X^\ := T x M n , x G M n ; 

(ii) and for every smooth function h : M n — > IR 

^£ / l* (0 COI*^ < TV{h) <c jr [ \x«\h)\ g dv g . 

An initial data u being given in L°°(M n ) fl BV(M. n ; dV g ), we want 
to find a solution u e = u e (t,x) to the initial value problem 

d t u e + dw g (f x (u e )) = eA/, x G M", t > 0, (4.1) 

u e (0,x) = u e (x), xGM™, (4.2) 
where A g denotes the Laplace operator on the manifold M n , 

A g v := div g gra.d g v 

n i d 2 v , dv v 
= q 3 ( r ■ ) . 



12 M. BEN-ARTZI AND P.G. LEFLOCH 

In (14.21) . Uq : M n — > M. is a sequence of smooth functions satisfying 

||Wq||lp(m«) < |ko||iP(M")) P e [Ij 00 ]) 

TV{ul) < TV(u ), 

e H u ollH 8 (M n ;dV ff ) < CTV(uo), for some C > depending only on M n , 
~^ u o a - e - on M n . 

(4.3) 

In the following, we use the notation u(t) instead of u(t, •). We begin 
with: 

Theorem 4.2. (Regularized conservation law.) Let f = f x (u) be a 
geometry- compatible flux on a Riemannian manifold (M ra ,g). Given 
u e Q G C°°(M n ) satisfying H4.3|) there exists a unique solution u e G 
C°°(R + x M") to the initial value problem fl4.1|) - fl4.2p . Moreover, for 
each 1 < p < oo the solution satisfies 

\\u e (t)\\ LP{M n. dVg) < \\u e (t')\\ LP{M n;dV g ), < t' < t (4.4) 

and, for any two solutions u e and v e , 

\\v'{t)-u e {t)\\ LHMn . dVg) < \\v e {t')-u'{t')\\ L i {Mn . dVg) , 0<t'<t. (4.5) 

In addition, for every convex entropy/ entropy flux pair (U, F x ) the so- 
lution u e satisfies the entropy inequality 

d t U(u e ) + div g (F(u e )) < e \U(u e ). (4.6) 

All proofs are postponed to the end of this section. 

Our next goal is to prove the strong convergence of the family {u e } e> o 
to an entropy solution of fl3.4l) - fl3.7p . as e — ► 0. We will assume here 
that the initial data has bounded total variation, and we refer to the 
forthcoming section for a more general argument which does not re- 
quire this assumption. From the estimates for the solutions obtained 
in Theorem 14.21 we can deduce a uniform total variation bound in space 
and time. 

Lemma 4.3. (BV bounds for diffusive approximations.) There exists a 
positive constant C\ depending on the geometry oflMJ 1 and ||mo||l°°(M' 1 )? 
only, such that the solutions u e given by Theorem \4-2\ satisfy 

TV(u e (t)) <e Clt (l + TV(u )), t G IR+, (4.7) 

\\dtu e (t)\\L 1 (M n -dv g ) < TV(uq) + e || A g ul\\L2(^ M n. dVg ^, t G R+. (4.8) 

Observe that the right-hand side of (14. 7p is independent of e, while 
the term e C\ \\ AgMgH^^n) arising in (14. 8p is dominated by CTV{uq) 
thanks to our assumption (14.31) . At this juncture, we refer to [1] in 
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which a sharp version of the above estimates is derived with a specific 
constant C\ depending on the geometry of the manifold. 

We are now in a position to state and prove the main result of this 
section. 

Theorem 4.4. (Convergence of the vanishing diffusion method on 
a Riemannian manifold.) Let f = f x (u) be a geometry-compatible 
flux on a Riemannian manifold (M n ,g). Given any initial condi- 
tion Uq G L°°(M n ) fl BV(M n ; dV g ) there exists an entropy solution 
u G L°°(]R + x M") to the initial value problem (I3.4p -( 13TT1) in the sense 
of Definition \3.S\ which is the limit of the sequence u e constructed in 
Theorem \4-2\ by vanishing diffusion. 
Moreover, this solution has the following properties: 

\\u{t)\\ L P(M";dV g ) < \\Uo\\ LP (M";dV g ), t G K+, 1 < p < OO, 

and, for some constant C\ > depending on ||wo||z,°°(M") an d the ge- 
ometry o/M n only, 

TV(u(t)) < e Cl * (1 + TV(u )), t G R+, 

\\u(t) - u(t')\\ L i (m n. dVg) < TV(u ) \t - t% 0<t'<t. 

Furthermore, if u, v are entropy solutions associated with some initial 
data u ,v , respectively, it hold 

\\v(t) - u(t)\\ L i {M n. dVg) < \\v - u \\ L i { M^;dv g ), t G R + . (4.10) 



Proof of Theorem By introducing a local chart of coordinates and 



using the condition f)3.5p . one can reduce the equation fl4.ll) to a para- 
bolic equation in the Euclidian space M": 

8 t tf + {d u fi){u*) d 3 u* = eg* (didju* - r* d k u<). (4.11) 

The local in time existence, uniqueness, and regularity of the solution 
u e follows as in the Euclidian setting. By putting together several 
charts to cover the whole manifold, one arrives at the local in time 
well-posedness for the diffusive conservation law (14.1 1) . Next, in order 
to extend the solution to arbitrary times, one needs a uniform estimate 
on the sup norm of the solution. This indeed follows in the form stated 
in (14. 4p for p = oo, from the standard maximum principle for parabolic 
equations. 

Next, in every strip [0, T] x M n the dual equation associated with 
(Oil reads 

dt¥ + 9x{f *x{u e ), grad ff <p) +eA g (p = 0, 
for which a "terminal value" problem with data given at t = T is 
now considered. This equation satisfies the maximum principle and, 
by duality, the inequality (I4.4p for p = 1 follows. The inequality for 
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intermediate values p G (l,oo) is then obtained by straightforward 
interpolation. 

We next proceed with the derivation of (14.51) . The function w e : = 
v e — u e satisfies the equation 

d t w e + div g (b e x w e ) = e A g w e , (4.12) 

where 

w Mv c )-f*(u e ) 
x ' v e - u e 

is a smooth vector field on M n . As in the argument above, the dual 
equation associated with (14.121) satisfies the maximum principle and 
so, by duality, 

lk e 0OIUi(M'w 9 ) < \\w e {t')\\ L i [m n. dVg)} <t' <t, 

which is precisely (14.51) . 

Finally, to show (14. 6 j) we multiply (14. ip by d u ll(u e ), where (U,F) is 
a convex entropy pair in the sense of Definition 13.31 We get 

d t U(u e ) + dw g F x (rf) = ed u U(u e )A g u e . 

In view of the definition of A g we can write in local coordinates 

d u U(u £ ) A g u e = A g U{u e ) - dlUiu') g jk d jU e d k u e . (4.13) 

Since d^U is non-negative (U being convex) and the matrix (gi ) is 
positive definite (since g is a Riemannian metric tensor), the second 
term in the right-hand side of (14.131) is non-positive. This establishes 
the desired inequality (14. 6p . □ 

Proof of Lemma \4-3[ To establish (14. 7p we rely on Lemma H~T1 without 
loss of generality, it is sufficient to estimate ||0i(^w e (t))||L 1 (M™;dv 9 ) f° r 
every smooth ip supported in a coordinate patch with ||'0||i, oo (M n ) < 1. 
Multiplying (14. lip by if) and differentiating with respect to x 1 , the 
function w e := di(if> u e ) is found to satisfy 

d t w e (t, x) + dj ((d u fl)(u f (t, x)) w e (t, x)) 
= A(x, u e , grad g u e )(t, x) + e A g w € (t, x), 

where A depends linearly on grad ff u e and smoothly upon x, u e (and 
is supported in 0). We now multiply the identity above by sgn(w e ) 
and integrate over M n . At this juncture we note that, since w e G 
H 2 (M n ;dV g ) at least, 

sgn(w e ) d t w e = d t \w t \, sgn(w e ) djiv e = dj\w e \, 
sgn(w e ) A g w e < A g |w e |, 
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where the latter inequality must be understood in the sense of distri- 
butions. 

Relying on the bound ||w e (i)||i,°°(M™) < ||^o||i oo (M n ) which follows 
from (14.41) . we then deduce that there exists a constant G\ > (de- 
pending on both ||mo||l°°(M") and the geometry of M") such that 

^ / [w'WldVgKQ + d [ \gmd g u%t)\ g dV g . 

This estimate can now be repeated for each of the (finitely many) vector 
fields . . . , X^ L \ and the desired conclusion (14. 7p finally follows 
from Gronwall's lemma. 

To establish (14.81) we differentiate (14. lip with respect to t and note 
that 

div 9 (d t f x (u e )) = d t divg (f x (u e )), 

A g (d t f x {u*)) = d t A g (f x (u e )). 

Also, the vector field <9i/ x (w e (t, x)) satisfies the condition (13.51) (with 
respect to the explicit dependence on x). 

It therefore follows that the function z e = d t {u e ) satisfies 

d t z%x) + JL((d u fi)(u%y))z%y)) y ^ = eA g z%x). (4.14) 

(In other words, the explicit dependence of (d u f^.)(u e (t,x)) on x is not 
differentiated.) Note that 

sgn(z^((d u fi)(u%y))z%y)) y=x 

= ^{{d u fi){u^x))\z^x)\) y=x . 

Defining the vector field r x := (d u f x )(u e (t,x)) \z e (t,x)\ we have, by the 
condition (13.51) . that the second term in the left-hand side of (14.141) 
satisfies 

d j(( d uf J x )(u £ (t,x)) \z e (t,x)\) = div g r x . 



Multiplying (14. 14p by sgn(2; e ) and integrating over M™, we obtain 

-/ 

dt jm" 

In turn, we can write 



\z £ {t)\dV g < 0. 



\z e (t)\dV g < / \z e (0)\dV g . 

M" JM n 
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Alternatively, this inequality can be derived from ( 14. 5 \} by choosing 
v(t,x) = v(t + a,x) with a — > 0. Using the equation (14.11) we can 
estimate the above term at t = 0: 

/ \z € (0)\dV g < ||grad fl Uol|Li(M»;dv 9 ) + e || A g Uo|| L 2 (M , Wg) 

< TV(u e ) + e || A 9 Wo|| L 2 (M n. dy!?) . 

Finally, applying this estimate to each of the vector fields . . . , X^ L \ 
we arrive at the estimate ( 14.81) . This completes the proof of Lemma fl~3l 

□ 



Proof of Theorem \4-4\ In view of the sup norm estimate (I4.4p and the 



uniform, space L 1 estimates ( 14.71) . we see that the sequence u e {t) is 
uniformly bounded in BV(M. n ; dV g ) for every time t and is therefore 
compact in L 1 (M n ; dV g ). Applying this argument at all rational times 
t and then picking up a diagonal sequence we can ensure that a subse- 
quence u ej (t) converges to some limit u(t) in the L l norm, as — > 0, for 
all rational t. Next, in view of the uniform time estimate (14.81) . the limit 
u(t) extends to all values of the time variable, with — > u(t) in 
L 1 (M n ; dV g ). Letting €j — > in the inequalities (14.61) one then deduces 
that u satisfies all of the entropy inequalities 

d t U(u) + dw g (F(u)) < 

in the weak sense ( I3.12p . 

In fact, by working in a localized coordinate patch and using the en- 
tropy formulation (14. 6 p with U{u) = \u — k\, F x {u) = sgn(u — k)(f x (u) — 
f x {k)), one can repeat the Kruzkov classical theorem in order to obtain 
a "localized" L 1 contraction property. By patching together (finitely 
many) coordinate patches we obtain a global estimate of the form 

\\v(t) -u(t)\\ L i {M n. dVg) <C\\v -u \\ L i ( M";dv g ), te[0,T], (4.15) 

where C > depends on T (and M n ). This is sufficient to imply unique- 
ness of the entropy solution u(t), hence the convergence of the whole 
family u e . At this point we can invoke the parabolic estimate ( 14.51) 
and conclude that C = 1, thus establishing the L 1 contraction prop- 
erty (14.101) . All the other estimates now follow from Theorem 14.21 and 
Lemma 14.31 □ 

5. The well-posedness theory 

Theorem assumes that the initial data has bounded variation and 
provides the existence of (locally BV) entropy solutions constructed by 
vanishing diffusion. In the present section, we obtain a generalization to 
L°° initial data and also establish the uniqueness of the entropy solution 
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in this larger class. We emphasize that the class of L°° solutions is 
completely natural for geometry-compatible conservation laws, since 
the L°° estimate is independent of the geometry while the BV estimate 
(14.91) . in general, depends upon it. 
Our generalization of Kruzkov's theory to manifolds is as follows. 

Theorem 5.1. (Well-posedness theory in L°° for geometry-compatible 
conservation laws.) Let f = f x (u) be a geometry- compatible flux on 
a compact, oriented, Riemannian manifold (M. n ,g). Given any ini- 
tial data u £ L 00 (M n ) there exists a unique entropy solution u £ 
L°°(]R + x M n ) to the initial value problem (13. 41) - 03.71) in the sense of 
Definition Iff. 31 Moreover, for each 1 < p < oo the solution satisfies 

\\ u {t)\\LP(M n ;dV g ) < H M o||LP(M";dV 9 ); t £ E + , (5.1) 

and, given any two entropy solutions u, v associated with some initial 
data Mo, vq, respectively, 

IK*) ~ u{t)\\ L i {m n. dVg) < \\v - WoIUhm";^), t £ K+. (5.2) 

Furthermore, we will see in the proof that the following inequality 
holds in the sense of distributions 

d t \v -u\ + dw g (sgn(u - v) (f x (v) - f x (u))) < 0. (5.3) 

Following DiPerna [3] we introduce the (larger) class of entropy 
measure- valued solutions and, in fact, establish a much stronger version 
of Theorem 15. II As observed in [3], Kruzkov's arguments take a simpler 
form in the measure-valued setting. We consider measure-valued maps 
v — v i,xi that is, weakly measurable mappings (t, x) £ M + x M n i— > v t x 
taking their values in the space of probability measures on K with 
support included in a fixed compact interval of R. The action of the 
measure v on a function U will be denoted by 

(v t , x ,U) := / U(u)du tjX (u). 

JM. 

The weak measurability property means that the map (yt,x-, is mea- 
surable in (t, x) for each U. 

Definition 5.2. Let f = f x (u) be a geometry- compatible flux on a Rie- 
mannian manifold (NL n ,g). Given any initial condition uq £ L°°(M' 1 ) ; 
a measure-valued map (t, x) £ M n x R + i— > i/ tjX is called an entropy 
measure- valued solution to the initial value problem (13. 4ft - (13. 71) if, for 
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every convex entropy /entropy flux pair (U, F x ) (see ( 13. lip ), 

[ [ ((vt,x, U) d t 9{t, x) + g x ((u t>x , F x ), grad g 9(t, x))) dV g (x)dt 



+ / U(u (x))9(0,x)dV g (x) > 0, 



M" 



(5.4) 

/or every smooth function 9 = 9(t, x) > compactly supported in 
[0,+oo) x M n . 

We will now prove that: 

Theorem 5.3. (The well-posedness theory in the measure- valued class 
for geometry-compatible conservation laws.) Let f = f x (u) be a geo- 
metry-compatible flux on a compact, oriented, Riemannian manifold 
(M n ,o). Let Mo be in L°°(M ri ) and v be an entropy measure-valued 
solution (in the sense of Definition 1 5.0) to the initial value problem 
(13 .41) - 03 .71) . Then, for almost every (t,x), the measure v t x is a Dirac 
mass, i. e. of the form 

Vt,x — 3u(t,x), 

where the function u G L°°(IR + x M n ) is the unique entropy solution 
to the problem (13. 4p - 03.71) in the sense of Definition \3.^\ Moreover, 



the solution satisfies the properties (15.11) - 05.31) . and the initial data is 
assumed in the strong sense 

limsup / \u(t, x) — Uq{x)\ dV g (x) — 0. (5.5) 
t-*o+ i M » 

We have already shown in Section 0] that, provided the initial data 
have bounded variation, entropy solutions can be constructed via van- 
ishing diffusion. The proof relies on the compactness the inclusion of 
BV into L l . Thanks to Theorem I5.3[ we can now provide an alter- 
native, more general proof which is valid for L°° initial data. In view 
of the uniform L°° estimate (i.e., the maximum principle which only 
requires the initial data to be in L°°), vanishing diffusion approxima- 
tions u e generate a measure-valued solution v. Moreover, v can be 
easily checked to satisfy all of the entropy inequalities and, in turn, 
Theorem 15.11 follows from Theorem 15.31 

Finally, we also have the following generalization of Theorem 14.41 

Corollary 5.4. The convergence result of the vanishing diffusion ap- 



proximations in Theorem \4-4\ remains valid if the initial data uq is 
solely in L°°(M n ). All estimates therein, except for the total variation 
estimate, still hold. 
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In the course of proving Theorem 15.31 we will need some a priori reg- 
ularity of measure-valued solutions, especially the fact that the initial 
data Mo is automatically assumed in a strong sense. 

Lemma 5.5. Let v = v ttX be an entropy measure-valued solution of 
( I3.4p - (13.7p . where it is assumed that m G L°°(M n ) and the supports of 
the measures v t ,x ore all included in a fixed common interval. Then: 

(a) For every convex function U = U(u) and every smooth function 
9 = 9(x), the function 



t^ / (u ttX ,U)6dV g (x) (5.6) 

has locally bounded total variation and admits a trace as t — > 0+. 

(b) For every function U = U(u, x), which is convex in u, measurable 
in x, and such that \U(u,x)\ < c\u\ + |t/(x)| with U G L 1 (M n ; dV g ) 
and c > 0, 

limsup/ (jy t ,x,U(.,x)) dV g (x) < / U (u (x) , x) dV g (x) . (5.7) 

(c) In particular, v assumes its initial data Uq in the following strong 
sense: 

limsup / / \u — uq{x)\ di*t,x{u) dV g {x) = 0. (5.8) 

Proof of Lemma 15.51 Using in the weak formulation (15 .4p a function 
9(t,x) = 9i(x)92{t), compactly supported in [0, oo) x M" and having 
9i,9 2 > 0, we obtain 



00 d9 



2 

dt 

oo 



(u, U) 6 X dV g {x)dt + 2 (O) / U(u ) 0i dV g (x) 

poo p 

>-/ 2 g x (gradg9 1 ,(is,F x ))dVg(x)dt 

JO JM n 

poo 

>-d / 9 2 dt, 



for some constant C\ > depending on 9\ (and the common support 
of Vt,x\ Thus the function 



V x {t) :=-Cit+ / {vt,x,U)9 l dV g 
satisfies the inequality 

poo m p 

~ Vi(t)-£dt<9 2 (0) U(u )9 1 dV g (x). (5.9) 
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Using a test-function 62 > compactly supported in (0, 00), we find 

That is, the function Vi(t) is decreasing and, therefore, has locally 
bounded total variation. Since it is uniformly bounded, Vi(t) has a 
limit as t — > 0+. This proves (a). 

To establish the item (b), we fix a time to > and consider the 
sequence of continuous functions 

9 e 2 (t)= l(t + e-t)/e, te[t Q ,to + e], 
(0, t>t Q + e. 

Relying on the property (a) above, we see that 

roc ine 

- / Vl (t)-ldt^V 1 (t +). 
Jo at 

Since #f(0) = 1 and t is arbitrary, (15. 9p yields 

VM = -Ct t + / (u^, U) 9t dV g {x) < [ U(u ) 6 1 dV 9 (x) 
for all t > and, in particular, for all Q\ = 6i(x) > 

lim / (u^tyOtdVgWK [ U{u )e x dV g {x). (5.10) 

*^°+ JM n Jm™ 

Note that the left-hand limit exists, in view of (a). 

Next, consider the set of all linear, convex and finite combinations 
of the form 

a, 61 j (x) Uj(u), 

j 

where a.j > 0, Ylj ocj — 1, the functions Uj are smooth and convex in u 
and the functions 0ij(x) > are smooth and compactly support, with 
moreover 

\Uj{u)9 ld {x)\ < c\u\ + \Uj(x)\ 

with c > and Uj G L 1 (M"; dV g ). This set is dense (for the uniform 
topology in u and the L l topology in x) in the set of all functions 
U = U(u, x) that are convex in u and measurable in x and satisfy 

\U(u,x)\ < c\u\ + \U(x)\ 

for some c > and U G L 1 (M n ; dV g ). Therefore, by density, (b) follows 
from (EH]). 

Finally, (c) follows from (b) by choosing U(u,x) = \u — uq(x)\. □ 
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Proof of Theorem \5.3[ In all of the following arguments, the inequali- 
ties should be understood in the sense of distributions. All steps can 
be justified rigorously by introducing test-functions in the usual way. 
Restricting attention in (15. 4p to functions 9 with compact support, we 
deduce that 

d t (u,U) +dw g (u,F) < (5.11) 
and so, after introducing the Kruzkov's entropies 

U(u,v) = \v- u\, F(u,v) = (f(v) - /(«)) sgn(v - u), 

we obtain 

d t (v,U(-,v)) + div 9 (u,F(-,v)) < 0, v eR. (5.12) 

Let fi be another entropy measure- valued solution to (13. 4p . We are 
going to combine (15.121) together with a similar statement for /i, that 
is, 

dt(n,U{u,-)) + dw g (fi,F{u, •)) < 0, ueR. (5.13) 
Introducing the tensor product v <g> jj, = v t)X ® \i t)X with, for instance, 

(vt,x ® Ht,xi U) := // U(u,v) dv tiX (u)d/i t , x (v), 

J JR2 

we can write (in the sense of distributions) 

d t (u ® fi, U) + div 9 (y ® /i, F) 

= (i/, ^(/i, f/) + div 9 (fjt, F)) + (fi, d t (u, U) + div 9 (u, F)) 

and we deduce from (I5.12p and (15.131) that, in the sense of distributions 

d t (v®fi,U) +div g (v® fx,F) < 0. (5.14) 

Next, integrating (15.141) over the manifold M n we find that, for all 
< f < t, 

/ (vt,x ® Vt,x,U) dVg(x) < / (v t ',x® Hv,x, U)dVg(x). (5.15) 

Letting t' — > in ( 15. 15ft and using that the two measure-valued solu- 
tions assume the same (Dirac-mass) initial data S Uo m at the time t = 
in the strong sense established in Lemma I5.5[ we find 

/ (vt< B ®th0,U)dV B (x)=O, t g R + . 

Therefore, {yt,x ® Ht.x,U} vanishes for almost every (t,x), which is 
equivalent to saying that v t x and n t)X coincide with the same Dirac 
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mass, say 6 u ( t ,x) f° r some value u(t,x). Since v and // are both en- 
tropy measure-valued solutions, the function u = u(t,x) is an entropy 
solution of the problem (I3.4l) - (l3.7p in the sense of Definition 13.31 

We establish some additional properties of the entropy solution u, as 
follows. Taking U(u) := \u\ p in (15.111) and integrating over the manifold 
M n yields 

- / U(u(t, x)) dV 9 (x) = - / (u ttX , U) dV 9 (x) < 0, 
at at Jm™ 

which leads to ( 15. ip . Using ( 15. 15ft with two distinct solutions v tjX = 
5 u (t,x) and fi t>x = 5 v (t,x) gives the L 1 contraction property (15.31) . This 
completes the proof of Theorem 15.31 □ 



6. General conservation laws and balance laws 

L 1 semi-group of entropy solutions on manifolds. Consider first 
the case of a conservation law posed on the one-dimensional torus T 1 = 
[0,1] 

d t u + \d x (kf(u)) — 0, u = u(t,x) e R, t G R+, x e [0,1], (6.1) 
k 

where k = k(x) is a given, positive function and / : R — > R is a convex 
function. In [TO] , a suitable generalization of Lax's explicit formula [7] 
was introduced for ( 16. lj) . When A; is not a constant, the characteristics 
of ( 16.11) are not straight lines but curves s i— >■ X(s) = X(s;y) (s > 0, 
i/Gl), given by 

d s X(s) = d u f(u(s,X(s))), X{0)=y. 

It was observed that along a characteristic the function v = v(s; y) : = 
u(s, X(s\ y)) is such that the "weighted flux" k(X) f{y) is constant. By 
introducing suitable left- and right- inverses of the function /, say f± x 
and then solving the equation k(X) f(v) = c, it follows that the whole 
family of all characteristic curves is described by 

d s X(s) = (d u f o ( ° ) , X(0) = y. (6.2) 

The following result was derived from a detailled analysis of this 
family of curves: 

Theorem 6.1. (See [8j [10].) The periodic, entropy solutions with 
bounded variation to the conservation law ( 16. ip are given by a gen- 
eralization of Lax's explicit formula via a minimization problem along 
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the curves (16.21) . Moreover, any two solutions u,v satisfy the L 1 con- 
traction property 

\\v(t) - u(t)\\ L i m < \\v(t') - u(t')\\ L i m , 0<t'<t. 

In view of the discussion in the previous sections, this result may 
seem surprising since the geometry-compatibility condition is not satis- 
fied here (except in the trivial case where k is a constant). Theorem 16. II 
motivates us to extend now our theory on Riemannian manifolds to 
general conservation laws that need not be geometry-compatible. 

From now on we consider a general conservation law (13 .4p associated 
with an arbitrary flux f x . First of all, we stress that the notion of 
entropy pair (U, F x ) should still be defined by the same conditions 
(13. lip as in the geometry-compatible case, but now we no longer have 
(I3.12j) . Instead, an entropy solution should be characterized by the 
entropy inequalities 

d t U(u) + div g (F(u)) - (div g F)(u) < (6.3) 

in the sense of distributions, for every entropy pair (U,F X ). The term 
(div 9 F) (u) is defined by applying the divergence operator to the vector 
field x — > F x (u), for every fixed u. 

Theorem 6.2. Let f = f x (u) be an arbitrary (not necessarily geometry- 
compatible) flux on a Riemannian manifold (Wl n ,g), satisfying the lin- 
ear growth condition (for some constant Co > 0) 

\f*(u)\ g < C {l + \u\), ueR,xeM n . (6.4) 

Then there exists a unique contractive, semi-group of entropy solutions 

uo G L 1 (M n ) i — ^ u(t) := S t u Q G L l {M n ) 

to the initial value problem (I3.4p - (l3.7p . 

The condition (16 .4p is required for the flux term (13 .4p to be an in- 
tegrable function on M n . Note that no uniform L p estimate is now 
available, not even in the L l norm, since the trivial function u = 
need not be a solution of the conservation law. The L 1 norm of a so- 
lution is finite for each time t, but generally grows as t increases. This 
result shows that the contraction property is more fundamental than all 
of the other stability properties derived earlier for geometry-compatible 
conservation laws. 

For general conservation laws, the stationary solutions u, determined 

by 

div g (f x (u(x))) =0, x G M", 
represent possible asymptotic states of time-dependent solutions. 
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Proof. The semi-group is constructed first over functions with bounded 
variation and then extended by density to the whole of L 1 . To this end, 
we re- visit the proof of Theorem I5.3[ and check that the key inequality 
(I5.14p . needed in the derivation of the L 1 contraction property, remains 
true regardless of the geometry-compatibility property. Note in passing 
that the stronger statement (15. lip is no longer valid under the present 
assumptions, and was precisely needed to ensure the stability in all LP 
spaces, which is no longer true here. 

It thus remains to establish (15. 14j) for solutions u, v of (13. 4p with 
bounded variation. We rely on standard regularity results (Federer 
[I], Volpert [12J) for such functions: pointwise values u±,v± can be 
defined almost everywhere with respect to the Hausdorff measure 3f n 
on R + x M n . These pointwise values coincide at points of approximate 
continuity. We compute the entropy dissipation measure 

/i := d t U(u,v) + div g (F(u,v)), 

and distinguish between the set Q U)V of points of approximate conti- 
nuity for both u and v, and the set § UtV of points where u has an 
approximate jump and v is approximately continuous, or vice- versa. 
The JCn-Hausdorff measure of the set N Ui „ := R + x M" \ G UjV U § u>v is 
zero. Using standard calculus for BV functions we can write 

V\eu,v ■= sgn(u - v) (d t u - d t v) + sgn(u - v) (div 9 (/(«)) - div 9 (f(v))) 
= 0. 

On the other hand, if B C $ u>v is a Borel subset consisting of points 
where (for instance) v is approximately continuous then 

fi{B) 

(n t \u+ -v\ + g x (n x ,sgn(u + - v) (f x (u+) - f x {v))\ 

- n t \u- - v\ - g x (n x ,sgn(u- - v) (f x {u-) - f x {v))jj dJ£ n (t,x), 

where n = (n t (t,x),n x (t,x) G R x T x M n is the propagation speed 
vector associated with the discontinuity in u. Now, by relying on the 
entropy inequalities, fiUt, x)} can easily be checked to be non-positive. 
In turn, the measure \i is also non-positive on $ UjV and, in turn, on the 
whole of R + x M n . The argument is complete. □ 

General balance laws. It is not difficult to generalize the well-posed- 
ness theory to the balance law 



d t u + div 9 (/(«, •)) = h(u, ■), 
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where / = f x (u,t) is a family of vector fields depending (smoothly) 
on the time variable t and on the parameter u, and h = h(u, t, x) is 
a smooth function. Our previous results should be modified to take 
into account the dependence of / in t and the effect of the source-term 
h. For instance, for L°° initial data, the L p estimate (15.11) should be 
replaced with 

||tt(t)||xP(M»;^) < C P (T) \\u{t')\\ LP{M n. dVg) + C'(T), < t' < t < T, 

where the constant C P (T) > 1 in general depends on T and C'(T) also 
depends upon g. Similarly, the contraction property should be replaced 
by 

\\v(t) - u(t)\\ L i(M»;dV g ) < C (T) \\v -U \\ L i {M n. dVg) , t E R+. 

for some constant Cq(T). We omit the details. 

7. Conservation laws on Lorentzian manifolds 

Globally hyperbolic Lorentzian manifold. Motivated by the math- 
ematical developments in general relativity, we now extend our theory 
to Lorentzian manifolds. 

Let (WL n+1 ,g) be a time-oriented, (n + l)-dimensional Lorentzian 
manifold, g being a pseudo-Riemannian metric tensor on M n+1 with 
signature (— , +,...,+). Recall that tangent vectors X on a Lorentzian 
manifold can be separated into time-like vectors (g(X,X) < 0), null 
vectors (g(X,X) = 0), and space-like vectors (g(X, X) > 0). A vector 
field is said to be time-like, null, or space-like if the corresponding 
property hold at every point. The null cone separates time-like vectors 
into future-oriented and past-oriented ones; it is assumed here that the 
manifold is time-oriented, i.e. a consistent orientation can be chosen 
throughout the manifold. 

Let V be the Levi-Cevita connection associated with the Lorentzian 
metric g so that, in particular, Vg = 0. The divergence div 9 operator is 
defined in a standard way which is formally similar to the Riemannian 
case. 

Definition 7.1. A flux on the manifold M™ +1 is a vector field x i— > 
f x (u) G T x WL n+1 , depending on a parameter u G R. The conservation 
law on (M. n+1 ,g) associated with f is 

div g (/(«)) =0, u: M n+1 -> R. (7.1) 

It is said to be geometry compatible if f satisfies the condition 

div g f x (u) = 0, uER, xe M n+1 . (7.2) 
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Furthermore, f is said to be a time-like flux if 

g x (d u f x (u),d u f x (u)) < 0, x G M n+ \ u G R. (7.3) 

Note that our terminology here differs from the one in the Riemann- 
ian case, where the conservative variable was singled out. 

We are interested in the initial- value problem associated with (17. ip . 
We fix a space-like hypersurface !Ko C M n+1 and a measurable and 
bounded function u defined on "Kq. Then, we search for a function 
u = u(x) G L°°(M n+1 ) satisfying (17. ip in the distributional sense and 
such that the (weak) trace of u on "Kq coincides with u : 

«[3<b = Uq. (7.4) 

It is natural to require that the vectors d u f x (u), which determine the 
propagation of waves in solutions of (17. ip . are time-like and future- 
oriented. 

We assume that the manifold M n+1 is globally hyperbolic, in the 
sense that there exists a foliation of M n+1 by space-like, compact, ori- 
ented hypersurfaces "K t {t G R): 

M n+1 = [J % t . 

Any hypersurface Ji to i s referred to as a Cauchy surface in M n+1 , while 
the family J{ t it G M) is called an admissible foliation associated with 
! K to - The future of the given hypersurface will be denoted by 

M+ +1 :=U^- 

Finally we denote by n* the future-oriented, normal vector field to each 
'Ht, and by g l the induced metric. Finally, along "K t , we denote by X 1 
the normal component of a vector field X, thus X 1 := g(X, n l ). 

Definition 7.2. A flux F = F x (u) is called a convex entropy flux asso- 
ciated with the conservation law (17. ip if there exists a convex function 
U : R -> R such that 

F x (u) = / d u U(u') d u f x {u') du', x G M n+1 , u G R. 
Jo 

A measurable and bounded function u = u(x) is called an entropy so- 
lution of the geometry- compatible conservation law (I7.ip - (l7.2p if the 
following entropy inequality 

/ g(F(u),grad g 8)dV g + [ g {F{u ), n ) 6^ dV go > 0. (7.5) 
Jm" +1 Juo 



GEOMETRY COMPATIBLE CONSERVATION LAWS ON MANIFOLDS 27 

for all convex entropy flux F = F x {u) and all smooth functions 9 > 
compactly supported in M™ +1 . 

In particular, in (17.51) the inequality 

div, (F(u)) < 
holds in the distributional sense. 

Theorem 7.3. Consider a geometry-compatible conservation law (17.11) - 
(17.21) posed on a globally hyperbolic Lorentzian manifold M n+1 . Let !Ko 
be a Cauchy surface in M n+1 , and Uq : "Hq — > W be a measurable and 
bounded function. Then, the initial-value problem (I7.ip - (l7.4p admits 
a unique entropy solution u = u(x) G L°°(M n+1 ). For every admis- 
sible foliation "Kt, the trace uji t exists and belong to L 1 (3-Ct), and the 
functions 

\\F (wj€ t |Ui(M t )) 

are non-increasing in time, for any convex entropy flux F. Moreover, 
given any two entropy solutions u, v, the function 

||/WJ-/'(<WIU W (7.6) 
is non-increasing in time. 

We emphasize that, in the Lorentzian case, no time-translation prop- 
erty is available in general, contrary to the Riemannian case. Hence, 
no time-regularity is implied by the L 1 contraction property. 

As the proof is very similar to the one in the Riemannian case, we 
will content with sketching the proof. Introduce a local chart 

x = (x a ) = (t, x j ), g := g a p dx a dx^, 

where by convention greek indices describe 0,1, ... ,n and latin indices 
describe 1, . . . ,n. By setting / =: (/"(«)) and using local coordinates, 
the conservation law (17. ip reads 

d a {\g x \ 1 / 2 f:(u(x)))=0, (7.7) 

where \g\ ■= det (gap)- Thanks to the assumption on /, for all smooth 
solutions (17.71) takes the equivalent form 

g x ((d u f x Ku(x)),gr a d g u(x)) := (d u f:){u(x)) (d a u){x) = 0. (7.8) 

In other words, setting x = (t,x) and / = (/*, /■?), 

(d u fl){u{t,x))d t u{t,x) + (d u fl)(u(t,x)) (d jU )(t,x) = 0, 

in which, since d u f is future-oriented and time-like, the coefficient in 
front of the time-derivative is positive 

d u f x (u) > 0, ueR,xe M n+1 . 
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To proceed with the construction of the entropy solutions, we add a 
vanishing diffusion term, as follows: 

dxv g (f(u c )) =eA/, (7.9) 

where A^ is the Laplace operator on the leaves K t of the foliation, that 
is in coordinates 

(d u f x )(u%x))d t u%x) + (d u fi)(u%x)) (d jU *)(t,x) 

= eg ij (t,x) (d&u 6 - F^d k u e )(t,x). 

In view of (I7.10p . we see that all of the estimates follow similarly as 
in the Riemannian case. For instance, multiplying ( 17.101) by a convex 
function U, the entropy inequality now takes the form 

d t (F°(u%x)) + d s {Fi{u e (t,Hf)) 

= d u U{u{t,x)) (dj°){u e {t,x))d t u*{t,x) 

+ d u U(u(t,x)) (d u fi)(u%x))(d jU e )(t,x) (7.11) 

= e d u U(u(t, x)) g ij (t, x) (d&u 6 - I* d k u € ) (t, x) 

= eAgU{u e ) -edlU{u £ ) g 3k d 3 u e d k u\ 

The metric g~ induced on the space-like leaves 'Kt is positive-definite 
and, therefore, the latter term above if non-positive and we conclude 
that, given a geometry compatible flux / and for every convex entropy 
flux F 

div g (F(u € )) = d t (F°(u € (t,x)) + dj{Fi{u e {t,x)) < eA 7 U(u e ). 

Note that the regularization (17.91) does depend on the specific folia- 
tion under consideration. However, by the contraction property (17.61) 
the limiting solution is unique and independent of the chosen regular- 
ization mechanism. 

Schwarzschild spacetime. In the context of general relativity, the 
Schwarschild metric represents a spherically symmetric empty space- 
time surrounding a black hole with mass m and is one of the most im- 
portant example of Lorentzian metrics. The outer communication re- 
gion of the Schwarzschild spacetime is a 1 + 3-Lorentzian manifold with 
boundary, described in the so-called Schwarschild coordinates (t, r, uj) 
by 

g = -(1 - — ) dt 2 + (1 - dr 2 + r 2 du 2 , 

with t > and r > 2m, while to describes the 2-sphere. There is 
an apparent (but not a physical) singularity in the metric coefficients 
at r = 2m, which corresponds to the horizon of the spacetime. This 
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spacetime is spherically symmetric, that is invariant under the group 
of rotations operating on the space-like 2-spheres given by keeping t 
and r constant. The part r 2 duo 2 of the metric is the canonical metric 
on the 2-spheres of symmetry. The spacetime under consideration is 
static, since the vector field d t is a time-like Killing vector. Moreover, 
this metric is asymptotic to the flat metric when r — > oo. Theorem 17.31 
extends to the exterior of the Schwarzschild spacetime, by observing 
that along the boundary r = 2m the characteristics of the hyperbolic 
equation are outgoing. 
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